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A  SURVEY  AND  ANALYSIS  OF 


Normal  Ionospheric  Absorption  Measurements 

OBTAINED  FROM  RADIO  PULSE  REFLECTIONS 


Larry  D.  Schultz  and  Roger  M.  Gallet 


ABSTRACT.  — An  analysis  of  published  midday  absorption  data  obtained  by  the  reflection 
method  of  vertical  incidence  sounding  quantitatively  establishes  the  existence  of  an 
equatorial  absorption  anomaly.  This  anomaly  is  the  abnormally  large  variation  of  the  ampli¬ 
tude  of  the  seasonal  variation  observed  at  equatorial  latitudes.  It  is  suggested  that  certain 
D-region  dynamical  processes  cause  this  anomaly.  In  addition,  the  analysis  disclosed  the 
unexpected  existence  of  a  significant  latitudinal  variation  of  absorption  not  due  to  the  cos  nx 
variation;  the  order  of  magnitude  of  this  latitudinal  variation  is  of  the  same  order  as  the 
variation  of  absorption  due  to  the  sunspot  cycle.  This  latitudinal  variation,  together  with 
the  latitudinal  variation  of  the  E-layer  critical  frequency,  indicates  that  the  electron  pro¬ 
ducing  mechanisms  for  the  lower  ionosphere  vary  systematically  with  latitude  and  solar 
activity.  For  middle  latitudes,  the  analysis  reconfirmed  results  obtained  previously  by  other 
authors. 

The  total  variation  of  lower  ionosphere  absorption  is  separated  into  the  four  categories: 
(1)  seasonal  variation  with  latitudinal  effects,  (2)  latitudinal  variation,  (3)  solar  cycle  varia¬ 
tion,  and  (4)  winter  anomaly  effects.  Empirically  derived  formulas  are  given  which  permit 
the  calculation  of  each  of  these  four  variations  and  also  the  calculation  of  the  total  lower 
ionosphere  absorption  losses  of  HF  radio  waves  propagating  via  the  ionosphere  both  at 
vertical  and  at  oblique  incidence. 


INTRODUCTION 

In  recent  years,  emphasis  in  the  literature  has 
been  on  the  effects  of  abnormal  ionospheric  absorp¬ 
tion  (i.e.,  polar-cap  and  auroral  absorption)  on  HF 
and  VHF  radio  waves  propagating  in  polar  regions. 
By  comparison,  the  effects  of  normal  ionospheric 
absorption  (i.e.,  nondeviative  D-region  absorption) 
on  propagation  have  been  somewhat  neglected,  even 
though  such  effects  apply  to  a  much  larger  area  of 
the  earth’s  surface.  This  neglect  of  normal  D-region 
absorption  effects  on  propagation  seems  due  to  the 
widespread  belief  that  all  absorption  effects  are 
well  known  and  well  understood,  and  that  little 
would  be  gained  by  additional  study  in  this  field. 


Such  a  belief  exists  because  different  independent 
analyses  of  absorption  data  obtained  at  middle 
latitudes  by  several  ionosphere  physicists  (Rawer 
1949  and  1952;  Piggott  1953;  Appleton  and  Piggott 
1954)  have  all  yielded  results  which  are  in  good 
agreement  with  each  other  and  with  results  pre¬ 
dicted  by  classical  photoionization  theory  (Appleton 
1937;  Nicolet  and  Bossy  1949;  Mitra  and  Jain  1963). 

The  usual  procedure  has  been  to  apply  the  classi¬ 
cal  theory  of  photoionization  to  equatorial  regions, 
because  photoionization  processes  are  assumed  to 
be  predominant  at  equatorial  regions  and  observa¬ 
tions  at  middle  latitudes  support  results  predicted 
by  the  theory.  This  procedure  is  not  correct,  how¬ 
ever,  since  results  predicted  by  the  theory  disagree 
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with  equatorial  observations.  The  reason  for  the 
disagreement  is  that  the  classical  theory  of  photo¬ 
ionization  incorrectly  assumes  that  the  D  region  is 
formed  by  uniform  photoionization  processes  and 
that  there  are  no  systematic  changes  with  latitude 
in  the  structure  of  the  D  region. 

An  early  indication  that  lower  ionosphere  proc¬ 
esses  vary  with  latitude  was  the  finding  of  a  syste¬ 
matic  latitudinal  variation  of  the  E-layer  critical 
frequency  ft:  in  addition  to  the  usual  variation 
caused  by  the  solar  zenith  angle.  This  variation  is 
given  by  (Harnischmacher  1950) 

fs  —  K  cos"  x> 

where 

K=  2.25+  1.5  cos  X  +  [0.01-0.007  cos  X]/?,3 
and 

=  0.21  +  0.12  cos  X+0.002 /?,*. 

The  quantity  \  is  the  solar  zenith  angle,  X  is  the 
geographic  latitude,  and  is  the  13-month  running 
average  of  the  mean  monthly  sunspot  number. 

Another  indication  that  D-region  absorption  proc¬ 
esses  were  not  the  same  for  all  latitudes  was  first 
reported  by  one  of  the  present  authors  (Delobeau 
and  Gallet  1954).  At  that  time,  it  was  pointed  out 
that  the  amplitude  of  the  seasonal  variation  of  ab¬ 
sorption  at  equatorial  latitudes  was  very  large  when 
compared  with  results  found  at  middle  latitudes 
and  also  with  results  predicted  by  classical  photo¬ 
ionization  theory.  It  was  suggested  that  this  anoma¬ 
lous  behavior  be  examined  further  when  more 
equatorial  absorption  data  became  available.  To 
the  best  of  the  present  authors’  knowledge,  such 
an  analysis  has  not  yet  been  given  in  the  literature. 

More  recent  evidence  that  D-region  formation 
processes  and,  therefore,  absorption  effects  vary 
with  latitude  has  come  from  the  study  of  galactic 
cosmic  rays  and  their  influence  on  the  ionosphere. 
These  phenomena  produce  large  amounts  of  elec¬ 
trons  (up  to  several  hundred  electrons/cm3)  at 
D-region  heights  in  high  geomagnetic  latitudes. 
Cosmic  rays  do  not  appear  to  have  a  diurnal  varia¬ 
tion  but  they  are  inversely  correlated  with  solar 
activity  (during  sunspot  maximum  the  sun  ejects 
more  plasmas  containing  magnetic  fields  which 
tend  to  repel  cosmic  ray  particles  entering  the 
solar  system).  The  effects  of  these  two  phenomena 
are  beyond  the  scope  of  this  paper,  and  therefore 
their  effects  are  not  discussed  or  included  in  this 
analysis. 


The  authors  became  interested  in  the  calculation 
of  lower  ionosphere  absorption  measurements 
through  the  problem  of  developing  a  computer 
method  for  calculating  HF  ionospheric  transmis¬ 
sion  losses.  In  the  course  of  the  development  of 
this  method,  a  search  of  the  literature  was  made  to 
find  suitable  formulas  that  would  permit  the  cal¬ 
culation  of  transmission  losses  due  to  lower  iono¬ 
sphere  absorption.  Since  the  primary  interest  was 
in  total  lower  ionosphere  absorption,  no  attempt  has 
been  made  here  to  separate  absorption  losses  into 
D-  and  E-region  components.  Such  a  separation  has 
been  attempted  by  others  (Bibl  and  Rawer  1951; 
Bibl,  Paul,  and  Rawer  1959;  Rawer  1960b;  Bibl, 
Paul,  and  Rawer  1961);  however,  theoretical  results 
from  these  separations  have  been  criticized  on 
experimental  grounds  (Fejer  1961).  The  authors 
feel  that  in  principle  there  are  certain  advantages 
in  separating  absorption  losses  into  various  compo¬ 
nents,  but  that  in  practice  the  separation  techniques 
are  not  yet  sufficiently  developed  to  permit  analyses 
similar  to  the  ones  performed  in  this  paper. 

The  search  of  the  literature  revealed  that  although 
there  has  been  much  work  on  the  theory  of  iono¬ 
spheric  attenuation  of  HF  transmissions,  little  has 
been  done  in  this  country  to  develop  formulas  and 
techniques  for  calculating  quantitatively  the  trans¬ 
mission  losses  involved.  Even  so,  of  the  few  formulas 
that  were  found  in  the  literature,  none  gave  results 
which  agreed  with  observations  made  at  equatorial 
latitudes.  Consequently,  to  produce  an  acceptable 
computer  method  for  calculating  transmission 
losses,  it  became  necessary  for  the  present  authors 
to  develop  absorption  formulas  valid  for  low  and 
middle  latitudes.  This  involved  the  analysis  of  a 
large  amount  of  original  published  absorption 
measurements  taken  at  stations  widely  separated 
in  latitude.  The  results  of  the  analysis  are  given  in 
this  paper. 

Finding  original  published  absorption  measure¬ 
ments  was  difficult  for  several  reasons.  First,  long 
series  of  continuous  absorption  measurements  taken 
over  one  or  more  solar  cycles  are  rare.  Second,  the 
majority  of  absorption  measurements  have  been 
made  at  middle  latitudes,  with  the  result  that  very 
little  data  exist  for  equatorial  latitudes.  In  addition, 
for  those  stations  for  which  long  periods  of  con¬ 
tinuous  data  existed,  the  data  were  not  always  avail¬ 
able  in  the  same  continuous  series  of  publications. 
In  the  case  of  Freiburg,  for  example,  absorption 
data  have  been  published  in  the  following  series: 
(1)  from  1948  to  1951  as  part  of  the  regular  monthly 
station  bulletin;  (2)  from  1952  to  1955  in  a  single 
report  (which  was  published  in  1956),  and  (3)  since 
the  year  1956  in  the  present  publication  form  of  a 
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separate  absorption  bulletin.  For  other  stations  the 
recording  of  absorption  measurements  has  not 
always  meant  that  the  observed  data  were  pub¬ 
lished.  In  the  case  of  Singapore,  for  example,  the 
authors  used  some  unpublished  data. 

Because  of  the  above  difficulties  encountered  in 
obtaining  the  published  values  of  the  original  meas¬ 
urements,  and  since  these  values  may  be  of  use  in 
future  studies,  all  of  the  data  analyzed  in  this  paper 
have  been  reproduced  in  the  appendix.  These  data 
are  the  published  monthly  median  midday  1  /Lvalues 
for  the  stations,  Slough,  Freiburg,  Singapore,  and 
Dakar.  The  latitudes  of  these  stations,  the  years  for 
which  data  were  published,  the  total  amount  of 
available  data,  and  the  range  of  sunspot  numbers 
(SSN)  over  which  data  were  observed  are  given 
below  in  table  1. 


Table  1.  —Equatorial  ionosphere  absorption  data 
used  in  this  report 


Station 

Latitude 

Period  of 
absorption 
measure¬ 
ments 

Years  of 
data 

available 

Range  of 
SSN 

(monthly 

mean) 

Slough 

51°31'N. 

1935-1953 

181 

0-201 

Freiburg 

48°06'N. 

1949-1961 

13 

0-254 

Dakar 

14°36'N. 

1951-1955 

41 

0-109 

Singapore 

1°19'N. 

1949-1954 

5 

0-144 

The  sources  of  the  above  data  are  given  in  the  ap¬ 
pendix,  page  29. 

The  purpose  of  this  paper  is:  (1)  to  analyze  and 
discuss  the  seasonal  and  solar  cycle  variation  of 
absorption  at  the  stations,  Slough,  Freiburg,  Dakar, 
and  Singapore;  (2)  to  examine  the  latitudinal  varia¬ 
tion  of  absorption  not  due  to  the  cos”  x  variation 
(X  is  the  solar  zenith  angle,  n  is  a  constant  for  a 
given  location);  (3)  to  develop  formulas  for  calculat¬ 
ing  HF  transmission  losses  due  to  lower  ionosphere 
absorption  which  separately  take  into  account 
seasonal  variation  of  absorption  with  latitudinal 
effects,  latitudinal  and  solar  cycle  variation  of  ab¬ 
sorption,  and  winter  anomaly  effects;  and  (4) 
to  give  in  one  publication,  for  the  convenience  of 
possible  future  studies,  all  of  the  midday  absorption 
data  analyzed  in  this  paper. 

In  this  paper,  a  brief  description  is  given  con¬ 
cerning  the  theory  of  lower  ionosphere  absorption 
and  the  procedure  for  taking  absorption  measure¬ 


1  In  the  case  of  Dakar,  all  of  the  ,4 -values  from  October  1951  through  December 
1955  were  obtained  at  1040  local  time. 


ments.  The  seasonal  variation  of  absorption  is 
discussed  in  detail  for  the  four  ionosphere  stations: 
Slough,  Freiburg,  Dakar,  and  Singapore.  Quanti¬ 
tative  formulas  are  given  for  calculating  the  seasonal 
variation  at  these  stations.  The  equatorial  and 
winter  anomalies  are  also  treated.  The  solar  cycle 
variation  of  absorption  is  given,  and  formulas  are 
derived  for  calculating  this  variation  as  a  function 
of  sunspot  number.  The  results  obtained  for  seasonal 
variation  of  absorption,  equatorial  and  winter 
anomalies,  and  the  solar  cycle  variation  are  com¬ 
bined  into  formulas.  These  results  are  compared 
with  published  results  obtained  by  other  authors. 
The  latitudinal  variation  of  absorption  not  due  to 
the  variation  in  cos  x  1S  discussed  next.  This  type 
of  variation  was  an  unexpected  result  of  the  data 
analysis,  and  to  the  best  of  our  knowledge,  has  not 
been  discussed  previously  in  the  literature.  An 
empirically  derived  formula  is  given  for  calculating 
transmission  losses  due  to  lower  ionosphere  ab¬ 
sorption  that  separately  takes  into  account  the 
effects  of:  (1)  seasonal  variation  with  latitudinal 
effects,  (2)  latitudinal  variation,  (3)  solar  cycle 
variation,  and  (4)  winter  anomaly  effects.  Finally, 
the  results  are  extended  to  oblique  and  nighttime 
absorption.  The  appendix  contains  the  published 
monthly  median  midday  values  of  the  absorption 
data  analyzed  in  this  paper,  together  with  the 
sources  of  the  data.  In  addition,  tables  of  various 
quantities  derived  during  the  course  of  the  data 
analysis  are  also  summarized  for  convenience  in  the 
appendix. 

LOWER  IONOSPHERE 
ABSORPTION  COEFFICIENT 
AND  ABSORPTION  INDEX 

Two  main  methods  are  generally  used  for  meas¬ 
uring  lower  D-  and  E-region  ionospheric  absorption. 
They  are  the  reflection  method  of  vertical  incidence 
sounding,  and  the  riometer  method  of  measuring 
changes  in  cosmic  noise.  This  paper  is  concerned 
only  with  absorption  measurements  obtained  from 
vertical  incidence  soundings. 

The  reflection  method  gives  the  apparent  reflec¬ 
tion  coefficient  p  of  the  ionosphere.  The  method  and 
the  operating  techniques  involved  are  discussed 
elsewhere  (Piggott  1953;  Piggott  et  al.  1957).  We 
now  give  a  brief  outline  of  the  classical  theory  for 
the  determination  of  the  D-region  absorption  by  the 
pulse  reflection  method.  The  Sen-Wyller  (1960) 
presentation  is  not  used  because  we  believe  that 
the  pulse  reflection  method  is  not  precise  enough 
to  warrant  this  more  sophisticated  treatment. 

The  real  refractive  index  pc  and  absorptive  index  x 
are  given  as  follows  (Budden  1961,  p.  172): 
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u2  v2=  1 

(1) 

2  = 

OJ  (jl)  b  V“ 

(2) 

where  (of)  is  the  angular  plasma  frequency,  a>  is 
the  angular  wave  frequency,  and  v  is  the  collision 
frequency.  Combining  equations  (1)  and  (2)  gives 
the  absorptive  index  x  as 


v  1  (^ 

*  2^(1 


-m2  +  x2)- 


(3) 


The  absorption  coefficient  k  occurs  in  the  exponent 
of  the  exponential  damping  factor  of  the  electric 
field  vector  of  the  radio  wave.  It  is  proportional  to 
the  absorption  index  and  is  given  by 


k=  (t o/c)x- 


(4) 


Using  the  approximation  (nearly  always  valid  in 
the  lower  ionosphere)  that  X“  —  M2  (Budden 
1961,  p.  173),  the  absorption  coefficient  is  given  by 
(Ratcliffe  1962,  p.  38)  as, 


v_  1~^ 
k  2c 


lire1  \  1 


Nv 


\mc€t)  /  pi  (a;  ±  to//) 2  - 


(5) 


(6) 


where  N  is  the  electron  density,  o)n  is  the  gyro- 
frequency^  is  the  velocity  of  light,  €o  is  the  dielectric 
constant  of  free  space,  and  e  and  m  are  the  electron 
charge  and  mass,  respectively.  The  plus  and  minus 
signs  in  equation  (6)  are  associated  with  the  ordi¬ 
nary  and  extraordinary  waves,  respectively. 
Throughout  the  remainder  of  this  paper,  we  con¬ 
sider  only  the  absorption  of  the  ordinary  wave. 

To  obtain  the  total  absorption  of  an  electromag¬ 
netic  wave  propagating  via  the  ionosphere,  it  is 
necessary  to  integrate  the  absorption  coefficient 
along  the  ray  path.  Thus,  for  a  radio  wave  emitted 
with  amplitude  E o  at  some  reference  distance  from 
a  transmitter,  the  amplitude  of  the  wave  reflected 
vertically  from  the  ionosphere  is  given  by 


Hi" 


k(z)dz 


(7) 


where  h'  and  h  are  the  virtual  and  true  heights  of 
reflection,  respectively  (Whitehead  1958;  Ratcliffe 


1962,  p.  402).  The  amplitude  £|  of  a  radio  wave 
reflected  from  the  ionosphere  is  also  given  by 
(Ratcliffe  1962,  p.  402) 


r  _ Eo 

E'~2h'p- 


(3) 


where  p  is  the  apparent  reflection  coefficient  of  the 
ionosphere.  Therefore,  the  amplitude  of  the  reflec¬ 
tion  coefficient  of  the  ionosphere  is  given  by 


p=  exp 


k(z)dz 


(9) 


(see  Budden  1961,  p.  325.  for  a  discussion  of  the 
phase  of  the  rtflec  on  coefficient).  In  units  of 
decibels,  the  attenuated  power  is  given  as  —20 
logio  p  which  is  equal  to  —8.7  log?  p. 

In  practice,  the  reflection  coefficient  is  measured 
each  day  at  noon  on  several  frequencies,  usually 
between  2  and  5  MHz.  The  method  and  the  operating 
techniques  involved  are  discussed  by  Piggott  (1953) 
and  Piggott  et  al.  (1957).  Since  the  absorption 
coefficient  is  proportional  to  Nv/co  +  ojh)2  if 
(a)±c oh)2>v2  (see  equation  (6)),  a  plot  of  logio  p 
versus  (co±c o//)2  will  show  whether  or  not  the  as¬ 
sumptions  underlying  equation  (8)  are  valid.  Employ¬ 
ing  these  considerations,  a  graphical  technique  is 
used  (Piggott  1953)  to  eliminate  those  values  of 
p  thought  to  be  influenced  by  the  occurrence  of 
deviative  absorption.  The  graphically  derived  re¬ 
sult  is  the  absorption  index  A,  which  is  independent 
of  frequency  and  which  gives  a  useful  index  of  the 
nondeviative  absorption  in  the  lower  ionosphere. 
The  absorption  index  A  is  related  to  the  reflection 
coefficient  p  by  the  relation 

A  (dB  MHz2)  —  8.7  (/+/w)2  loge  p,  (10) 

where  fn  is  the  gyrofrequency  (Appleton  1937). 

The  absorption  index  A  is  used  to  obtain  the  ab¬ 
sorption  at  vertical  incidence  on  other  frequencies 
not  too  near  the  critical  frequency  of  an  ionized 
region  (i.e.,  for  frequencies  at  which  deviative 
absorption  is  not  significant).  The  formula  used  to 
calculate  this  transmission  loss  La  in  decibels 
is  given  by 


La(dB)=AI(f+fHy 


(ID 


where  A  is  the  absorption  index. 

For  oblique  incidence  propagation,  the  value  of 
absorption  obtained  from  equation  (11)  is  multi- 
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plied  by  a  factor  secant  </>,  where  0  is  the  angle  of 
incidence  at  the  midpoint  of  the  absorbing  region. 
In  using  this  simple  procedure,  however,  caution 
should  be  exercised  in  situations  involving  multiple 
hops,  because  the  absorption  index  is  usually  ob¬ 
tained  only  from  first  order  reflections  from  the 
E  and  F  regions.  Later  in  this  paper,  the  modi¬ 
fication  of  the  vertical  incidence  formula  for  the 
case  of  oblique  incidence  is  discussed  in  greater 
detail. 

To  calculate  the  absorption  for  times  other  than 
noon,  equation  (11)  is  multiplied  by  a  diurnal 
variation  factor  (cos  0.893x/cos  0.893xi2)rf.  The 
quantity  X12  is  the  solar  zenith  angle  in  degrees  at 
local  noon,  and  x  is  the  solar  zenith  angle  at  the 
time  of  the  calculation.  The  factor  0.893  determines 
the  sunset  and  sunrise  times  at  a  height  of  100  km 
above  the  earth.  The  experimental  value  of  the  ex¬ 
ponent  d  is  of  the  order  of  unity  (Taylor  1948)  with 
values  ranging  from  0.6  to  1.3  being  reported 
(Appleton  and  Piggott  1954).  The  reported  values, 
however,  are  based  upon  a  sunrise-sunset  factor 
of  unity  rather  than  0.893. 

The  general  expression  for  the  absorption  losses 
of  an  obliquely  incident  wave  propagating  via  the 
ionosphere  is  given  by 


L„(dB)  =  v4(dB  MHz'2) 


cos  0.893x  ~\d  secant 
.cos  0.893X12J  ”(/+//, )2’ 


(12) 


where  the  symbols  have  their  previously  defined 
meanings.  All  of  the  quantities  in  equation  (12) 
except  A ,  can  be  calculated  from  the  geometry  of 
the  ray  path.  The  index  A  depends  upon  the  electron 
density  and  collision  frequency  profiles,  and  there¬ 
fore,  A  must  be  determined  empirically.  This 
determination  is  made  in  the  next  section.  In  the 
latter  part  of  this  paper,  the  denominator  term  con¬ 
taining  the  frequency  is  modified  to  take  into 
account  absorption  losses  when  the  angular  wave 
frequency  co  is  comparable  to  the  collisional  fre¬ 
quency  v . 


THE  GEOGRAPHIC  AND 
TEMPORAL  VARIATION  OF 
THE  ABSORPTION  INDEX 

The  measured  midday  values  of  the  absorption 
index  A  exhibit  a  short-term,  annual  variation  and 
a  long-term,  11-year  variation.  The  short-term  varia¬ 
tion  is  due  to  the  seasonal  change  of  the  solar  zenith 
angle,  and  the  long-term  variation  is  due  to  the 
sunspot  cycle.  These  two  variations  are  illustrated 


in  figures  1  through  4.  In  these  figures,  the  published 
monthly  median  midday  2  values  oi  A  are  shown  for 
the  stations,  Slough,  Freiburg,  Singapore,  and 
Dakar.  In  addition,  the  13-month  smoothed  /Lvalues 
and  the  13-month  smoothed  sunspot  numbers  are 
also  given;  these  are  designated  Al3  and  i?13,  re¬ 
spectively.  For  convenience  the  /Lvalues,  and  in 
addition  the  A  ^-values  are  given  in  the  appendix. 

The  annual  variation  of  the  /Lvalues  (i.e.,  varia¬ 
tion  of  the  nondeviative  absorption)  is  due  to  the 
seasonal  changes  in  the  solar  zenith  angle.  For 
Slough  and  Freiburg,  maximum  absorption  occurs 
during  the  summer  months  (figs.  1  and  2),  but  for 
Singapore  and  Dakar,  maximum  absorption  occurs 
slightly  after  the  overhead  transit  of  the  sun.  Be¬ 
cause  these  stations  are  at  equatorial  latitudes, 
maximum  absorption  occurs  twice  a  year  (figs.  3 
and  4). 

The  11 -year  variation  of  absorption  is  related  to 
the  11-year  variation  of  the  sunspot  cycle.  This  solar 
cycle  variation  is  illustrated  in  figures  1  through  4 
by  the  R  13-curves.  The  year-to-year  variation  of  the 
/Lvalues  is  obtained  as  a  running  average  of  the 
^-values  over  a  seasonal  cycle  (13  months).  This 
year-to-year  variation  is  also  shown  in  figures  1  to 
4.  It  can  be  seen  that  the  correlation  between  the 
curves  of  absorption  and  solar  activity  is  very  good. 

Thus,  we  see  that  absorption  values  for  each 
ionosphere  station  undergo  two  different  variations: 
(1)  a  seasonal  variation  and  (2)  a  solar  cycle  varia¬ 
tion.  These  variations  are  now  determined 
quantitatively. 

ANALYSIS  OF  THE  SEASONAL 
VARIATION  OF  ABSORPTION 

Observations  indicate  (Best  and  Ratcliffe  1938; 
Rawer  1949)  and  theory  suggests  (Appleton  1937; 
Nicolet  and  Bossy  1949)  that  the  seasonal  variation 
of  absorption  is  given  by 

A  =A13K  cos"  x  (13) 

where  A 13  is  the  13-month  running  average  of  the 
^-values,  x  is  the  solar  zenith  angle,  and  K  and  n 
are  constants  for  a  given  station. 

The  parameters  K  and  n  in  equation  (13)  are 
normally  determined  in  the  following  manner.  From 
the  published  A -values  one  obtains  first  the  13- 
month  running  average  of  the  A -values  (this  is  the 
average  over  one  period  of  the  seasonal  variation). 

-*  In  the  case  of  Dakar,  the  published  lvalues  from  October  1951  through  December 
1955  were  obtained  at  1040  local  time.  Noon  values  were  estimated  by  multiplying  the 
published  values  by  the  diurnal  variation  factor  cos  x/cosx„M„  =  1.065,  taking  the  ex¬ 
ponent  d=  1.0  in  this  case. 
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Normal  Ionospheric  Absorption  Measurements 


FIGURE  3.  — Variation  of  the  published  monthly  median  midday  values  of  the  absorption  index  A.  The  quantity  A1S  is  the  smoothed, 
13-month  average  of  the  ^-values.  /?13  is  the  smoothed  13-month  average  of  the  sunspot  numbers.  The  occurrence  of  the  double 
maxima  usually  corresponds  to  the  semiannual  zenith  passages  of  the  sun. 


Next  each  A -value  is  divided  by  its  corresponding 
,4 13-value;  that  is,  for  example,  the  January  1950 
,4 -value  is  divided  by  the  v4i3-value  for  January 
1950.  The  ratio  A\A  i3  is  designated  /.  Thus,  a  set  of 
/-values  is  obtained,  one  for  each  month.  Finally, 
the  average  of  all  the  /-values  for  each  month  is 
calculated.  These  average  /-values  are  designated 
7.  The  quantities  A ,  A i3,  /,  and  I  are  given  in  tables 
A-I  through  A-IV  of  the  appendix.  For  convenience 
the  /-values  are  also  summarized  below  in  table  2. 

The  variations  of  I  with  cos  x  f°r  the  different 
stations  are  shown  in  figures  5  through  8.  For  the 
middle  latitude  stations,  Slough  and  Freiburg, 
one  notices  that  a  Large  variation  in  cos  x  results  in 
a  relatively  small  variation  in  absorption.  More 
precisely,  for  Slough  a  240  percent  change  in  cos  x 
results  in  a  63  percent  change  in  7,  and  for  Freiburg 
a  184  percent  variation  in  cos  x  corresponds  to  a 
43  percent  change  in  7.  Thus,  at  middle  latitudes, 
the  variation  of  absorption  (i.e.,  the  /-values)  is 
of  the  order  of  one-fourth  the  variation  of  cos  x- 

In  view  of  the  above  results,  one  would  expect  a 
small  amplitude  for  the  seasonal  variation  of  ab¬ 


Table  2 .  —  Summarization  of  I-  values 


I— A!  A  ,;1 

Slough 

Freiburg 

Singapore 

Dakar 

Jan. 

0.990 

0.876 

0.940 

0.822 

Feb. 

0.886 

0.901 

0.973 

0.699 

Mar. 

0.881 

0.881 

1.052 

0.806 

Apr. 

1.138 

1.086 

1.089 

1.225 

May 

1.208 

1.142 

1.023 

1.326 

June 

1.180 

1.140 

0.889 

0.770 

July 

1.211 

1.160 

0.865 

1.067 

Aug. 

1.147 

1.104 

1.004 

0.954 

Sept. 

1.035 

1.063 

1.085 

1.299 

Oct. 

0.826 

0.838 

1.140 

1.197 

Nov. 

0.740 

0.806 

1.035 

0.838 

Dec. 

0.884 

0.876 

0.910 

0.896 

sorption  at  equatorial  latitudes,  since  the  seasonal 
variation  of  cos  x  is  small  at  the  Equator  (i.e.,  the 
midday  sun  is  almost  always  near  the  zenith).  How¬ 
ever,  instead  of  the  small  amplitude  variation  of  the 
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FIGURE  4.  — Variation  of  the  published  monthly  median  midday  values  of  the  absorption  index  A.  The  quantity  A 13  is  the  smoothed, 
13-month  average  of  the  /T values.  R 13  is  the  smoothed  13-month  average  of  the  sunspot  numbers.  The  occurrence  of  the  double 
maxima  usually  corresponds  to  the  semiannual  zenith  passages  of  the  sun. 


order  of  one-fourth  that  of  the  variation  of  cos  x 
observed  at  middle  latitudes,  one  finds  an  amplitude 
variation  of  approximately  3V2  times  the  cos  x  varia¬ 
tion.  For  Dakar,  for  example,  a  change  in  cos  x  °f 
28  percent  is  accompanied  by  a  change  of  89  per¬ 
cent  in  /;  and  for  Singapore  the  change  in  cos  x  is 
only  9  percent,  whereas  the  variation  of  I  is  of  the 
order  of  32  percent. 

The  abnormally  large  variation  of  equatorial 
absorption  is  somewhat  unexpected,  since  the 
midday  solar  radiation  is  almost  constant  at  these 
latitudes.  This  anomalous  behavior  was  pointed  out 
by  one  of  the  present  authors  several  years  ago 
(Delobeau  and  Gallet  1954).  At  that  time,  however, 
data  were  insufficient  to  determine  quantita¬ 
tively  the  magnitude  of  the  variation  (i.e.,  the  value 
of  the  exponent  n  in  the  seasonal  factor  cos”  x)- 
With  the  data  presently  available,  the  calculation 
of  the  magnitude  of  the  equatorial  seasonal  variation 
was  performed.  The  results  of  these  calculations 
are  discussed  later  in  the  text. 

In  addition  to  the  double  maxima  and  minima 
seen  in  figures  7  and  8  for  Singapore  and  Dakar, 


one  also  notices  in  each  graph  that  both  absorption 
maxima  and  one  of  the  minima  lag  the  cos  x  curves 
by  approximately  18  days  (as  determined  from  a 
cross-correlation  analysis).  This  lag  is  of  the  order 
of  the  magnitude  of  the  thermal  lag  observed  just 
below  D-region  heights  (Reed  1962),  which  is  in  all 
probability  associated  with  the  dynamics  of  the 
D  region.3  To  take  this  delay  into  account  for 
determining  the  seasonal  variation  of  absorption, 
it  was  necessary  to  shift  the  values  of  /  in  figures 
7  and  8  by  an  amount  corresponding  to  18  days. 
The  values  resulting  from  this  translation  are  the 
ones  used  in  the  subsequent  analysis. 

To  determine  the  seasonal  variation  of  absorption 
for  each  station,  standard  logarithmic  plotting  pro¬ 
cedure  is  used  to  transform  the  graphs  of  figures 
5,  6,  7,  and  8  into  the  form  shown  in  figures  9  through 
12.  In  these  latter  graphs,  the  systematic  seasonal 
changes  of  absorption  are  readily  seen  by  following 
the  arrows  from  one  month  to  the  next  (the  numbers 

3  D-region  dynamics  also  account  for  the  abnormally  large  variations  in  the  amplitude 
of  the  equatorial  absorption  mentioned  above.  A  discussion  of  the  mechanisms  involved 
in  D-region  dynamics  is  given  later. 


10 


Normal  Ionospheric  Absorption  Measurements 


Figure  6.  — Variation  of  absorption  with  cos  For  a  middle- 
latitude  station,  a  large  variation  in  cos  \  results  in  a  small 
variation  of  absorption. 


Figure  5. —  Variation  of  absorption  with  cos  x •  For  a  middle- 
latitude  station,  a  large  variation  in  cos  x  results  in  a  small 
variation  of  absorption. 


in  the  figures  represent  months  of  the  year— 1  for 
January,  2  for  February,  etc.).  Figures  9  through  12 
were  plotted  as  I  versus  V2  cos  \  in  order  to  center 
the  data  in  the  graph. 

The  straight  lines  in  figures  9  through  12  were 
determined  by  least-squares  curve-fitting  tech¬ 
niques.  In  the  case  of  Slough  and  Freiburg,  only  the 
nonwinter  months  of  March  through  October  were 
used  in  determining  the  equations  of  the  lines  be¬ 
cause  of  the  winter  anomaly  phenomena  discussed 
below.  Using  the  Student  t  coefficient  and  assuming 
that  the  data  are  independent  and  normally  dis¬ 
tributed,  the  values  of  K  and  n  were  calculated  (see 
equation  (13))  for  the  four  stations  and  are  given 
below  in  table  3  together  with  95  percent  confidence 
limits  for  n.  The  lines  in  figures  9  through  12  were 
extrapolated  beyond  cos  x  equal  to  unity  in  order  to 
show  the  large  changes  in  slope  of  the  seasonal 
variation  for  the  different  stations. 


Figure  7.  — Variation  of  absorption  with  cos  x-  For  an  equatorial 
station,  a  small  variation  in  cos  x  results  in  a  large  variation  of 
absorption  (compare  wjth  figs.  5  and  6).  Notice  the  pronounced 
lag  of  the  maxima  of  I  relative  to  the  maxima  of  cos  x-  There 
is  also  a  slight  lag  in  the  minimum  at  July.  The  overall  lag  is  of 
the  order  of  18  days  (determined  by  a  cross-correlation  analysis). 
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Figure  8.— Variation  of  absorption  with  cos  For  an  equatorial 
station,  a  small  variation  in  cos  \  results  in  a  large  variation 
of  absorption  (compare  with  figs.  5  and  6).  Although  not  as 
apparent  as  at  Singapore  (fig.  7),  an  overall  lag  of  approximately 
18  days  also  occurs  for  this  station. 


Table  3 .  —  Empirical  values  of  constants  K  and  n 
I  =  K  cos  tt  x 


Station 

K 

n 

95%  confidence 
limits  for  n 

Slough . 

1.35 

0.72 

0.57  ^  n  ^  .88 

Freiburg . 

1.25 

0.66 

0.46  ^  n  ^  .85 

Dakar . 

1.10 

1.53 

0.64  2.43 

Singapore . 

1.10 

2.40 

1.43  ^  ng  3.37 

Classical  photoionization  theory  predicts  a  value 
of  1.5  for  the  exponent  n  for  a  simple  Chapman- 
type  absorbing  region  (Appleton  1937).  For  more 
realistic  models  of  the  D  region,  values  of  n  less 
than  1.5  are  always  obtained  (Nicolet  and  Bossy 
1949;  Appleton  and  Piggott  1954;  Mitra  and  Jain 
1963).  Therefore,  the  values  of  n  given  in  table  3 
for  Slough  and  Freiburg  are  not  in  disagreement 
with  the  theory.  However,  for  the  equatorial 
stations  Dakar  and  Singapore,  the  observations 
are  in  disagreement  with  the  theory— particularly 
in  the  case  of  Singapore. 


Probably  the  reason  that  classical  photoionization 
theory  does  not  explain  the  abnormal  amplitude  of 
the  seasonal  variation  of  absorption  at  the  Equator 
is  because  the  theory  does  not  take  into  account 
dynamic  processes  of  the  D  region  caused  by  geo¬ 
graphical  factors,  meteorological  effects,  systematic 
seasonal  variations  of  the  ionosphere  for  different 
latitudes,  and  similar  effects.  These  processes 
give  rise  to  large-scale  motions  of  air  masses 
between  the  northern  and  southern  hemispheres. 
These  motions  generally  follow  with  a  lag,  the  lati¬ 
tudinal  position  of  the  sun  (Hess  1959).  This  hori¬ 
zontal  motion  produces  turbulence  which  results 
in  a  mixing  of  the  atmospheric  constituents  in  the 
vicinity  of  the  Equator.  As  a  result  of  this  mixing, 
electrons  in  the  lower  ionosphere  are  brought 
down  to  lower  heights.  This  mixing  results  in  an 
increase  in  the  amount  of  absorption,  since  absorp¬ 
tion  is  proportional  to  /  vNedh,  where  v  is  the  colli¬ 
sion  frequency  and  N v  the  electron  density;  and 
since  v  increases  exponentially  for  decreasing 
heights  in  the  ionosphere.  (The  above  description 
is  necessarily  qualitative  for  two  reasons.  First, 
quantitative  knowledge  of  atmospheric  turbulence 
at  D-region  heights  and  the  factors  producing  it 
are  not  available.  Second,  there  is  a  scarcity  of 
equatorial  absorption  measurements  with  which 
to  compare  observations  with  theory.  Consequently, 
effects  of  turbulence  on  absorption  are  beyond  the 
scope  of  the  present  paper.) 

For  Slough  and  Freiburg  (figs.  9  and  10),  the 
winter  months  of  November,  December,  January, 
and  February  give  proportionately  more  absorption 
than  the  nonwinter  months.  This  is  the  well-known 
‘"winter  anomaly”  phenomena  (Appleton  and 
Piggott  1954).  To  account  for  the  “winter  anomaly” 
when  using  equation  (13),  it  is  necessary  to  use  an 
empirical  winter  effect  coefficient  W .  This  coefficient 
is  unity  for  all  months  except  November,  December, 
January,  and  February  (the  months  of  the  “winter 
anomaly”).  The  value  of  W  for  these  months  is 
such  that  the  observed  values  of  absorption  are 
obtained  instead  of  the  values  given  by  the  straight 


Table  4.  —  Average  winter  effect  coefficient ,  W,  for 
Slough  and  Freiburg 


Month 

Slough 

Freiburg 

Average 

November . 

1.14 

1.18 

1.16 

December . 

1.64 

1.48 

1.56 

January . . . 

1.69 

1.39 

1.54 

February . 

1.18 

1.17 

1.18 
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Figure  9.— Seasonal  variation  of  absorption  as  a  function  of  cos  X-  The  nonwinter  months  of  March  through  October  are  seen  to  lie 
very  close  to  the  straight  line  (for  these  months  the  correlation  coefficient  equals  0.80).  The  winter  months  deviate  markedly  from 
the  straight  line;  this  is  due  to  the  winter  anomaly  effect. 


lines  in  figures  9  and  10.  The  JF-values  of  these  four 
months  are  summarized  in  table  4  and  are  also 
included  in  figures  9  and  10. 

Introducing  the  winter  effect  coefficient  W  results 
in  the  following  expression  for  the  seasonal  variation 
of  the  A -values: 


A—AizKW  cosn  x  (14) 

where  for  each  of  the  stations  the  values  of  K  and  n 
are  given  in  table  3,  and  for  the  middle  latitude 
stations,  the  values  of  W  are  given  in  table  4. 
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Figure  10,  — Seasonal  variation  of  absorption  as  a  function  of  cos  X-  The  nonwinter  months  of  March  through  October  are  seen  to  lie 
very  close  to  the  straight  line  (for  these  months  the  correlation  coefficient  equals  0.96).  The  winter  months  deviate  markedly  from 
the  straight  line;  this  is  due  to  the  winter  anomaly  effect. 
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Figure  11.  — Seasonal  variation  of  absorption  at  Dakar. 
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FIGURE  12.  — Seasonal  variation  of  absorption  at  Singapore. 
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ANALYSIS  OF  THE  SOLAR  CYCLE 
VARIATION  OF  ABSORPTION 

In  this  section,  the  analysis  of  the  variation  of 
absorption  with  the  sunspot  cycle  is  discussed. 
This  variation  of  absorption  with  the  solar  cycle 
is  shown  by  plotting  the  smoothed  absorption  index 
values  A\w  as  the  ordinate  with  the  smoothed  sun¬ 
spot  values  7?  13  as  the  abscissa.  The  resulting  graphs 
are  shown  in  figures  13  through  16.  The  Slough 
absorption  series  covered  two  solar  cycles  and  was 
plotted  for  convenience  on  two  graphs  (figs.  13 
and  14). 

For  Slough  and  Freiburg  (figs.  13,  14,  and  15),  one 
notices  that  the  absorption  values  are  different  for 
the  ascending  and  descending  phases  of  the  solar 


cycle.  In  the  case  of  Slough,  the  magnitude  of  the 
values  are  also  seen  to  change  from  one  solar  cycle 
to  the  next  (compare  figs.  13  and  14).  A  similar 
behavior  has  been  observed  in  the  relationship 
between  the  critical  frequency  in  the  F>  region, 
foF2,  and  sunspot  number  (Ostrow  and  PoKempner 
1952). 

To  simplify  as  much  as  possible  the  final  expres¬ 
sion  for  the  absorption  index  A ,  it  was  decided  to 
use  a  single  linear  relation  for  all  of  the  Slough 
data  and  a  single  linear  relation  from  all  of  the 
Freiburg  data.  In  the  case  of  Singapore  (fig.  16), 
one  straight  line  is  seen  to  give  a  good  fit  to  all  of 
the  observed  data.  However,  for  Dakar  (fig.  16), 
there  are  insufficient  data  to  determine  the  nature 
of  the  relation  between  absorption  and  sunspot 


VARPATION  OF  ABSORPTION  WITH  SUNSPOT  NUMBER 
 -SLOUGH  [CYCLE  17)- 


ALL  DATA  N935- 19531 


V+44 


1*0.00502  R, 


1  +  0.00fl04  R 


FicURK  13.  — Absorption  as  a  function  of  solar  activity  for  sunspot  cycle  17.  The  straight  lines  give  the  linear  relation  of  absorption 
with  sunspot  number.  The  dashed  line  gives  the  relation  for  cycle  17.  The  solid  line  gives  the  relation  for  the  combined  data  of 
cycles  17  and  18  (see  fig.  14). 
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Figure-  14.  —  Absorption  as  a  function  of  solar  activity  for  sunspot  cycle  18.  The  straight  lines  give  the  linear  relation  of  absorption  with 
sunspot  number.  The  dashed  line  gives  the  relation  for  cycle  18.  The  solid  line  gives  the  relation  for  the  combined  data  of  cycles 
17  and  18  (see  fig.  13). 


number.  For  Dakar,  only  the  value  Ao  of  absorption 
at  sunspot  minimum  can  be  deduced  — and  then 
only  approximately. 

The  expression  relating  absorption  and  sunspot 
number  is  assumed  to  be  linear  and  is  given  by 

A13=Ao[1  +  bRn]  (15) 

The  values  of  A{)  and  b  for  the  four  ionosphere  sta¬ 
tions  analyzed  are  given  below  in  table  5  and  are 
also  given  in  the  appropriate  graphs  of  figures  13 
through  16. 


Table  5 .  —  Values  of  coefficients  A0  and  b 


Station 

Ao 

b 

Slough . 

326 

0.00502 

Freiburg . 

304 

0.00550 

Dakar  . 

270 

Singapore . 

289 

0.00443 
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Figure  15.  — Absorption  as  a  function  of  solar  activity.  Notice  the  difference  in  absorption  for  the  ascending  and  descending  phases  of 
the  solar  cycle. 
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Figure  16. —  Absorption  at  Dakar  and  Singapore  as  a  function  of  solar  activity. 
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SUMMARY  OF  STATION  ANALYSES 

In  summarizing  the  results  of  the  analyses  for 
four  ionosphere  stations,  the  combined  seasonal 
and  solar  cycle  variation  of  the  median  midday 
/I -values  is  given  by  (see  equations  (14)  and  (15)) 

A  =  A»K[l  +  bRl*]W  cos”*,  (16) 

where  the  numerical  quantities  are  given  in  table  6. 

The  winter  effect  coefficient  W  is  summarized  in 
table  7. 


TABLE  6.  — Empirical  equations  for  ionospheric 
absorption  at  four  stations 


Station 

Latitude 

A  ~A{)K  [1  +  bR  I*]  W  oos"x 

Slouch 

51.5°N. 

A  =  440  [1+5.02  cos0*72* 

Freiburg 

48.1°N. 

A  =  380  [1  +  5.50  (10)  :1 rcosll-'i8x 

Dakar 

14.6°N. 

A  =  295  [1  +  .  .  .  R  ,.i]  W  cos1-*3  x 

Singapore 

i 

1.3°N. 

A  =  318  [1+4.43  (10)-;,«,3]r  c.»s2Wx 

COMPARISON  WITH  PUBLISHED  RESULTS 
OBTAINED  AT  OTHER  STATIONS  BY  OTHER 
AUTHORS 


In  this  section  the  results  of  analyses  by  other 
authors  are  presented.  These  results  were  obtained 
from  the  literature  on  D-region  absorption  studies. 
Unfortunately,  the  amount  of  material  dealing  with 
seasonal  variation  is  small,  probably  because  most 


Table  1. —  Winter  effect  coefficient ,  W 


For  the  months 

Slough 

Freiburg 

Dakar 

Singapore 

March  through 
October . 

1  00 

1.00 

1.00 

1.00 

November . 

1.14 

1.18 

1.00 

1.00 

December . 

1.64 

1.48 

1.00 

1.00 

January . 

1.69 

1.39 

1.00 

1.00 

February . 

1.17 

1.17 

1.00 

1.00 

D-region  studies  have  been  concerned  with  the 
diurnal  variation  of  absorption  (not  the  concern  of 
this  paper).  Since  the  present  authors  did  not  have 
access  to  the  other  authors’  original  measurements, 
only  the  results  of  their  statistical  analyses  can  be 
given. 

The  results  of  these  analyses  are  summarized 
below  in  table  8. 

The  results  in  table  8  are  in  general  agreement 
with  the  present  authors  results  which  are  sum¬ 
marized  in  table  6. 


LATITUDINAL  VARIATION  OF 
ABSORPTION 

The  coefficient  A{)K  and  exponent  n  appearing  in 
formula  (16)  should  be  independent  of  latitude  if: 


Table  S.  — Analyses  of  D-region  absorption  by  other  authors 


Station 

I  atitude 

n 

bx  103 

Period 

References 

Ibadan 

7.4  °N. 

0.90  (Sunspot  MAX) 
2.00  (Sunspot  MIN) 

2.6  (/=2.4  Mc/s) 

3.4  (/=  5.7  Mc/s) 

1953  to  1958 

Skinner  and  Wright  1961 

Tsumeb 

26.1  °S. 

1.5 

1957  to  1958 

Umlauft  1961 

Delhi 

28.6  °N. 

0.80  for  D  and  E 
region 

0.62  for  D  region 

1958  to  1959 

Ran,  Mazumdar,  and  Mitra  1962 

Prince  Rupert 

54.3  °N 

0.50 

1949  to  1958 

Davies  and  Hagg  1955 

Churchill 

58.8  °N. 

0.30 

1955  to  1966 

Peebles  1956 

Baker  Lake 

64.3  °N. 

0.10 

1955  to  1956 

Peebles  1956 
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(1)  the  D  region  is  formed  by  uniform  photoioniza¬ 
tion  processes;  and  (2)  the  structure  of  the  at¬ 
mosphere  at  D-region  heights  does  not  change 
systematically  with  latitude.  Examination  of  table 
6  shows,  however,  that  both  the  coefficient  A0K  and 
the  exponent  n  vary  systematically  with  latitude. 

The  variation  of  A{)K  with  latitude  is  illustrated 
in  figure  17.  Assuming  that  the  variation  is  linear 
(he.,  essentially  assuming  that  Slough  and  Freiburg 
correspond  to  one  observation,  and  that  Singapore 
and  Dakar  correspond  to  another  observation),  the 
relation  between  A{)K  and  latitude  k  is 

Ai)K= 286[1  -f-  0.0087X]  (17) 

where  k  is  the  geographic  latitude  expressed  in 
degrees. 


The  variation  of  n  with  latitude  is  illustrated  in 
figures  18  and  19.  Because  of  the  limited  amount  of 
data  available  at  equatorial  latitudes,  the  confidence 
limits  for  n  are  very  large  for  Singapore  and  Dakar. 
By  comparison,  the  confidence  limits  for  Slough 
and  Freiburg  are  very  narrow.  The  nature  of  the 
variation  of  n  with  latitude  cannot  be  given  with 
good  accuracy,  primarily  because  of  the  large  con¬ 
fidence  limits  for  the  equatorial  stations.  As  can  be 
seen  from  figures  18  and  19,  both  a  linear  relation 
and  a  cosine  relation  fit  the  observed  values  equally 
well.  The  results  illustrated  in  these  figures  indi¬ 
cate  the  inadvisability  of  using  a  single  value  of  n 
for  all  latitudes  — as  is  the  practice  of  most  predic¬ 
tion  services  (Rawer  1960a).  We  have  chosen  to  use 
the  linear  relation  because  of  its  simplicity,  and 
because  it  gives  a  cutoff  in  the  seasonal  variation 


ATITUOINAL  VARIATION  OF  THE  QUANTITY  AQK 


SLOUGH 


^FREIBURG 

haiui.Mui-M-  ilj 


; .  • _ : _ 
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t  IOURE  17.—  Latitudinal  variation  of  the  amplitude  of  absorption  excluding  the  variation  due  to  cos  x- 
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beyond  the  Arctic  Circle.  The  linear  variation  of  n 
is  driven  by 

=  2.25  — 0.032X  (18) 

where  X  is  the  latitude  expressed  in  degrees. 

The  linear  and  cosine  curves  given  in  figures  18 
and  19  were  obtained  by  a  least-squares  fit  of  the 
data  from  Singapore,  Dakar,  Freiburg,  and  Slough. 


In  addition  to  the  values  for  these  stations,  the 
values  obtained  by  other  authors’  analyses  (see 
preceding  section)  are  also  displayed  in  figures  18 
and  19.  As  can  be  seen,  these  other  results  are  not 
in  disagreement  with  the  present  authors’  results; 
and  with  the  exception  of  the  values  for  Ibadan  and 
Delhi  at  sunspot  maximum,  all  of  the  values  are  in 
good  agreement  with  each  other. 


LATITUDINAL  VARIATION  OF  n  IN  I  ^KCOSnX 


SINGAPORE 


[BADAN 
SS  MINI 


DAKAR 


TSUMEB 


IBADAN 


SLOUGH 


DELHI 


PRINCE  RUPERT 


CHURCHILL 


FREIBURG 


SAKER  LAKE 


CONFIDENCE  LIMITS 
{95%) 


40  50 

X  (LATITUDE) 


FIGURE  18.  — Variation  of  the  exponent  n  in  the  seasonal  factor  cos"x  as  determined  from  the  values  from  Singapore,  Dakar,  Freiburg, 
and  Slough.  The  solid  line  gives  the  assumed  linear  variation  with  latitude.  The  dashed  lines  give  the  locus  of  values  of  two  standard 
deviations  from  the  solid  line  (i.e.,  ±2S«|a)-  The  intervals  on  each  of  the  station  points  are  95  percent  confidence  limits. 
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VARIATION  OF  n  IN  WKCOSnX 


SINGAPORE 


i  EBAD4N 
HfSS  MTN) 


. .  !  i . ! 


SUUlEBt- 


[MKflRj 


1:<  -fog  cos  \ 

■  A — . - 1  — * — 


[IBADAN 

l(SSMAX) 


- - A - 1 

•4$  LOUGht- 


PRINCE  RUPERT 


DELHI 


CHURCHILL 


FREIBURG 


iCONROENCE  LIMITS  FT 


(95%) 


BAKER 


Figure  19. —  Variation  of  the  exponent  n  in  the  seasonal  factor  cos"  x  as  determined  from  the  values  from  Singapore,  Dakar,  Freiburg, 
and  Slough.  The  solid  line  gives  the  assumed  cosine  variation  with  latitude.  The  dashed  lines  give  the  locus  of  values  of  two  standard 
deviations  from  the  solid  line  (i.e.,  ±2S«]x).  The  intervals  on  each  of  the  station  points  are  95  percent  confidence  limits. 


At  the  present  time,  the  latitudinal  effects  of  the 
winter  anomaly  are  not  known.  The  JF-factors  dif¬ 
ferent  from  unity  appearing  in  tables  6  and  7  can 
probably  be  applied  only  within  a  certain  range  of 
geographic  latitudes.  In  using  the  JF-values,  the 
following  procedure  is  recommended:  (1)  apply  the 
values  only  to  latitudes  above  30°,  (2)  use  the  average 
of  the  Slough  and  Freiburg  values  for  latitude  60°, 


and  (3)  let  these  average  values  decline  linearly 
north  and  south  of  60°  such  that  they  have  the  value 
unity  at  latitudes  30°  and  90°.  Use  a  corresponding 
procedure  during  local  winter  in  the  southern 
hemisphere. 

There  does  not  seem  to  be  any  systematic  varia¬ 
tion  with  latitude  of  the  coefficient  b  (see  table  5). 
Consequently,  the  use  of  a  constant  value  of  6,  say 
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6  =  0.005,  for  all  locations  i<*  recommended,  or  the 
appropriate  station  value,  whichever  is  more 
convenient. 

Table  9. —  The  winter  anomaly  factor ,  W,  for 
stations  at  or  north  of  30°  latitude  for  all  months 


Month 

Latitude  (X) 

W 

Nov. 

30  ^  X  ^  60 

1. +  0.0083  (X  — 30) 

Dec. 

30  ^  X  ^  60 

1. +0.0282  (X  — 30) 

Jan. 

30  X  60 

1. +0.0269  (X  — 30) 

Feb. 

30  ^  X  ^  60 

1. +  0.0089  (X  — 30) 

Nov. 

60  X  90 

1. +0.0083  (90  — X) 

Dec. 

60  ^  X  ^  90 

1. +0.0282  (90  — X) 

Jan. 

60  ^  X  ^  90 

1. +0.0269  (90 -X) 

Feb. 

60  ^  X  ^  90 

1. +  0.0089  (90 -X) 

All  other  cases . 

1.00 

EMPIRICALLY  DERIVED  FORMULA  FOR 
CALCULATING  TRANSMISSION  LOSSES 

In  deriving  an  overall  formula  for  the  calculation 
of  the  absorption  index  A  which  takes  seasonal, 
solar,  and  latitudinal  variations  into  account, 
a  certain  amount  of  accuracy  is  lost  as  compared 
to  using  individual  formulas  for  each  station. 
In  spite  of  this  disadvantage,  however,  there  are 
numerous  situations  which  warrant  the  use  of  a 
general  formula.  With  this  limitation  in  mind,  we 
suggest  the  following  formula  for  the  total  median 
lower  ionosphere  absorption  at  midday  of  a  HF 
ionospheric  propagation 

LaU\\>,)=  (19) 

where  /  and  fi  are  given  in  MHz  and  where  A 
is  given  by 


A( dB)  =  286  [1^+0. 0087X]  [  1  +  0.005^]  W  cos^.  (20) 

latitudinal  solar  cycle  winter  seasonal 
variation  variation  anomaly  variation 

effect  with 

latitudinal 

effect 


where  \  and  X12  are  in  degrees,  n  —  2.25  —  0.032X 
and  W  is  given  in  table  9. 

To  calculate  absorption  for  times  other  than 
midday,  equation  (10)  is  multiplied  by  a  diurnal 
variation  factor  [cos  0.893x/cos  0.893xi2]d.  This 
factor  is  discussed  on  page  5  of  the  text. 

Equation  (20)  differs  from  other  formulas  found 
in  the  literature  in  the  following  points:  (a)  there 
is  an  explicit  latitudinal  variation  which  does  not 
depend  on  cos  x  and  (b)  the  exponent  of  the  seasonal 
variation  factor  cos  x  is  a  function  of  A  (rather  than 
a  constant  for  all  latitudes). 


CALCULATION  OF  OBLIQUE 
INCIDENCE  ABSORPTION 

Modification  of  Vertical  Incidence  Formula 

For  long-distance  telecommunications,  it  is 
necessary  to  calculate  the  absorption  of  radio  waves 
propagating  obliquely  via  the  ionosphere.  As  stated 
on  page  4,  it  is  the  usual  practice  to  multiply  the 


vertical  incidence  absorption  by  a  secant  factor 
that  is  proportional  to  the  angle  of  incidence  of 
the  wave  at  the  absorbing  region.  This  procedure 
would  seem  justified  if  the  radio  wave  penetrated 
entirely  the  absorbing  region;  however,  this  situa¬ 
tion  does  not  always  obtain  at  lower  HF.  Another 
difficulty  that  occurs  at  lower  HF  is  that  the  wave 
frequency  and  the  collision  frequency  are  com¬ 
parable.  The  net  result  is  that  oblique  incidence 
absorption  measurements  below  about  4  MHz  do 
not  follow  an  inverse  square  frequency  law  variation; 
that  is,  absorption  proportional  to  the  inverse  second 
power  of  frequency.  At  the  present  time,  it  is  not 
known  how  much  of  the  observed  discrepancy  is 
due  to  the  assumption  of  a  secant  of  the  angle  of 
incidence  variation  and  how  much  is  due  to  the 
assumption  of  an  inverse  frequency  variation. 
Since  the  collision  frequency  is  not  known  to  within 
a  factor  of  about  two,  one  normally  uses  an  effective 
collision  frequency.  The  present  justification  for 
the  use  of  an  effective  collision  frequency  is  that  it 
compensates  empirically  for  most  of  the  discrepan¬ 
cies  between  observed  and  predicted  oblique 
absorption  values. 
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The  previous  considerations  suggest  that  the 
vertical  incidence  formula  (19)  should  be  modified 
for  oblique  incidence  absorption  so  that  the  formula 
becomes: 


r  /  m\  y4(dBMHz2)  sec 

La{dB)-  v*tw- +{f+fHy 2 


(21) 


where  v  is  the  effective  collision  frequency.  The 
additive  term  u2/4772  in  the  denominator  of  equation 
(21)  is  significant  when  v 2  is  approximately  equal  to 
(co-f  a)#)2,  which  occurs  at  lower  HF.  Values  of  the 
effective  collision  frequency  obtained  from  oblique 
incidence  data  range  from  18.3  X 106  sec-1  to 
20.1  X  (10)6  sec-1  (Josephson  1966;  Lucas  and  Hay- 
don  1966).  The  present  authors  determined  a  theo¬ 
retical  value  of  20.2  X  (10)6  sec*1  as  follows.  A  plane 
electromagnetic  wave  is  assumed  incident  upon  the 
ionosphere  given  by  the  electron  density  profile 
and  the  collision  frequency  profile  shown  in  figure 
20.  The  assumed  electron  density  model  holds  for 


COLLISION  FREQUENCY,  (sec-1) 

Figure  20.  — Electron  density  and  collision  frequency  profile 
used  for  the  calculation  of  the  absorptive  index. 


approximately  sunspot  maximum,  and  the  collision 
frequency  distribution  is  a  fifth  degree  polynomial 
approximation  to  Nicolet’s  1963  values.  Using  the 
Appleton-Hartree  equation,  the  absorption  index 
(see  page  4)  is  calculated  as  a  function  of  height 
for  a  frequency  of  3  MHz  and  is  plotted  in  figure  21. 
The  absorption  index  has  its  maximum  value  at  a 
height  of  63  kilometers,  and  this  corresponds  to  a 
collision  frequency  of  20.2  X  (10)6  sec-1  (i.e.,  effec¬ 
tive  collisional  frequency).  For  a  plane  wave  at  30 
MHz,  the  absorptive  index  peaks  at  58  kilometers, 
and  this  corresponds  to  an  effective  collisional 
frequency  of  40X(10)6  sec*1  (fig.  21).  At  a  wave 


frequency  of  30  MHz,  however,  the  effect  of  the 
additive  term  u2/4tt2  is  negligible  (less  than  5  per¬ 
cent  change  in  the  calculated  absorption  values). 


Figure  21.  — Calculated  absorptive  index  as  a  function  of 
frequency  and  height. 


Comparison  Between  Calculated  and  Observed 
Oblique  Absorption  Values 

Equations  (20)  and  (21)  were  used  to  calculate 
absorption  values  over  the  path  between  Long 
Branch,  Illinois,  and  Boulder,  Colorado,  at  midday 
for  the  months  November  1958  through  November 
1959.  Two  sets  of  values  were  obtained  correspond¬ 
ing  to  effective  collision  frequency  values  of  0  and 
20  X  (10)6  sec-1,  and  these  values  are  plotted  in 
figure  22.  The  curve  labeled  A  corresponds  to  a 
collisional  frequency  of  zero  (i.e.,  use  of  equation 
(19),  and  the  curve  labeled  B  was  obtained  from 
equation  (21)  with  a  collisional  frequency  of 
20  X  (10)6  sec-1.  During  the  observational  period, 
an  oblique  HF  and  VHF  propagation  experiment  was 
conducted  over  this  path  (Blair,  Davis,  and  Kirby 
1961),  and  included  in  the  experiment  was  a  trans¬ 
mission  loss  study  at  5  MHz,  using  a  radiated  power 
of  39  kw.  Open  circuit  voltage  was  recorded  at  the 
receiver  and  the  hourly  median  decibel  values  of 
voltage  above  one  microvolt  were  obtained.  The 
system  loss  excluding  the  ionospheric  loss  (Barg- 
hausen  et  al.  1969)  was  calculated  and  compared 
with  an  analysis  of  the  5  MHz  data  (Davis  1969)  to 
obtain  the  measured  monthly  median  midday  ab¬ 
sorption  values.  These  absorption  measurements 
are  also  plotted  in  figure  22  (solid  curve).  As  can 
be  seen,  the  agreement  between  the  calculated 
absorption  and  the  observed  absorption  appears  to 
be  satisfactory.  Unfortunately,  observations  at  5 
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MHz  were  not  available  throughout  the  entire  period 
of  the  experiment  because  of  equipment  difficulties. 

As  we  have  indicated  in  this  paper,  our  formula 
was  derived  from  vertical  incidence  data  and  then 
was  generalized  from  theoretical  considerations 
to  the  calculation  of  oblique  data.  In  contrast  to 
this  approach,  a  formula  derived  from  considerations 
of  only  oblique-incidence  absorption  data  is  given 
by  Lucas  and  Haydon  (1966), 

j  ,  ,m_677(l  +  0.0037fl13)  (cos  0.881x)13  sec  4>  _ 

L(dbj"  10.2  +(/+/«)*•»  *  (22) 

We  have  plotted  as  curve  C  in  figure  22,  for  compari¬ 
son  purposes,  the  absorption  values  for  the  Long 
Branch-Boulder  path  as  calculated  by  equation 
(22).  As  can  be  seen,  this  latter  formula  on  this 
particular  circuit  does  not  agree  with  the  observa¬ 
tions  as  well  as  results  calculated  by  equations 
(20)  and  (21).  Other  comparisons  between  observed 
and  calculated  signal  strengths  also  indicate  satis¬ 
factory  agreement  when  equations  (20)  and  (21) 
are  used  (Barghausen  et  ah  1969).  The  diurnal  varia¬ 
tion  factor  previously  discussed  has  the  solar  zenith 
X  multiplied  by  the  factor  0.893  to  take  into  account 
layer  sunrise  and  sunset  times  at  100  km.  However, 
in  the  study  indicated  in  figure  22,  this  factor  was 
changed  to  0.881  when  the  results  of  calculations 
by  equations  (20),  (21),  and  (22)  were  compared. 


Figure  22.  — Comparison  of  calculated  and  observed  5-MHz 
absorption  on  the  Long  Branch-Boulder  path  (1292  km). 


Nighttime  Absorption 

Nighttime  absorption  at  upper  HF  is  not  notice¬ 
able,  but  for  oblique  propagation  at  lower  HF  it 
amounts  to  a  few  decibels.  A  study  of  nighttime 
field  strengths  (Barghausen  et  al.  1969)  indicated 
that  the  quantity 


(1.0  + 0.0005/?13)  cos'- X12  [-0?,nQn8Q!3x)J  W 

[_cos  (0.893xi2)J 

appearing  in  equation  (20)  (multiplied  by  the 
diurnal  variation  term  discussed  on  page  5), 
when  set  to  the  value  0.01,  gives  calculated  results 
in  agreement  with  observations.  In  the  above  ex¬ 
pression,  X12  represents  the  solar  zenith  angle  at 
midday,  and  the  factor  0.893  takes  into  account 
sunrise  and  sunset  times  at  a  height  of  100  km  above 
the  earth. 


CONCLUSION 

The  seasonal  and  solar  cycle  variation  of  iono¬ 
spheric  absorption  is  analyzed  at  each  of  the 
stations,  Slough,  Freiburg,  Singapore,  and  Dakar. 
For  the  middle  latitude  stations,  Slough  and  Frei¬ 
burg,  the  analysis  reconfirms  the  results  of  other 
published  analyses,  and,  in  addition,  gives  quantita¬ 
tively  the  effects  of  the  winter  anomaly  for  the 
northern  winter  months.  In  the  case  of  the  equatorial 
stations,  anomalous  absorption  is  found  to  occur. 

A  comparison  of  the  results  of  the  four  stations 
reveals  the  unexpected  existence  of  a  significant 
latitudinal  variation  of  absorption  not  due  to  the 
cos"  x  latitudinal  variation  (x  is  the  solar  zenith 
angle;  n  is  a  constant  for  a  given  location).  The 
latitudinal  variation  is  of  the  same  order  of  magni¬ 
tude  as  the  solar  cycle  variation  of  absorption. 

The  total  variation  of  lower  ionosphere  absorption 
is  separated  into  the  four  categories:  (a)  seasonal 
variation  with  latitudinal  effects,  (b)  latitudinal 
variation,  (c)  solar  cycle  variation,  and  (d)  winter 
anomaly  effects.  Empirically  derived  equations  are 
given  which  permit  the  calculations  of  each  of  these 
four  variations  and  also  the  calculation  of  the 
transmission  losses  of  HF  radio  waves  propagating 
via  the  ionosphere. 
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Appendix 


Absorption  Data  Sources 

This  section  gives  the  sources  of  the  published 
original  absorption  data  obtained  by  the  reflection 


Absorption  data  sources 


Station 

Period 

Publication 

Slough 

Jan.  1935- June  1953 

Monthly  Bulletin  of  Ion¬ 
osphere  Characteris¬ 
tics,  B-Series, 

Slough,  England 

Freiburg 

Dec.  1948-Mar.  1951 

Jan.  1952-Dec.  1955 

Jan.  1956- June  1962 

Observations  Iono- 
spheriques,  SPIM-0 
Series 

Observations  de 

F  Absorption  Iono- 
spherique,  SPIM-R16 
Ionospharen  Daten, 
Absorptions  Mes- 
sungen,  Nr.  1 
through  Nr.  7 

Singa¬ 

pore 

Sept.  1949-June  1953 

July  1953-Apr.  1954 

Monthly  Bulletin  of  Iono¬ 
sphere  Characteris¬ 
tics,  B-Series, 

Slough,  England 

Private  Information 

Dakar 

May  1951- Dec.  1955 

Observations  de 
l’Absorption  Iono- 
spherique,  SPIM- 
R16 

method.  For  the  convenience  of  future  studies,  and 
in  view  of  the  difficulty  of  obtaining  the  absorption 
material  which  was  dispersed  in  a  large  number  of 
different  bulletins,  it  is  believed  valuable  to  assemble 
in  homogeneous  form  all  the  monthly  median  A- 
values  from  the  beginning  of  the  series  — January 
1935  until  the  last  published  value,  June  1962.  The 
sources  of  this  data  are  given  in  the  table. 

With  the  exception  of  part  of  the  Singapore  data, 
all  of  the  data  listed  in  the  above  table  and  used 
in  this  report  was  made  available  to  the  present 
authors  through  the  facilities  of  the  Aeronomy  and 
Space  Data  Center  of  the  Environmental  Science 
Services  Administration  in  Boulder,  Colorado. 


Description  of  Tables 

The  appendix  contains  four  tables,  one  table  for 
each  of  the  ionosphere  stations,  Slough,  Freiburg, 
Dakar,  and  Singapore.  Each  table  contains  the 
following: 

1.  The  published  absorption  index  ^-values 

from  the  sources  given  in  the  preceding 
table. 

2.  The  smoothed  absorption  index  ,4 13-values. 

3.  The  ratio  I  =  A/A  13  for  each  month  for  which 

data  were  available. 

4.  The  average  /-values  for  each  month, 

designated  /.  These  values  are  also  given 
in  table  2  of  the  text. 
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Appendix 


Table  A-L  —  Slough 

YEAR 

lg 35  1  1936  1  1937  1  1938  |  1939  I  1940  i  19M  '  I  1942  I  1943  I  1944  | 

A -VALUES 


JAN, 

194 

_ 

400 

380 

410 

- 

824 

342 

^298 

356 

JAN. 

FEB. 

3  6b 

378 

434 

3ff- 

34o 

484 

388 

322 

234 

294 

FEB. 

MAR. 

370 

282 

470 

364 

446 

L30 

384 

380 - 

348 - 

356 

MAR. 

APR. 

446 

438 

1022 

540 

658 

512 

408 

492 

416 — 

550 

APR, 

MAY 

610 

868 

826 

564 

634 

598 

388 

388 

- 468 — 

526 

MAY 

JUNE 

_ 

.  636 

934 

534 

564 

635 

4^6 - 

3§l - 

436" 

544 

JUNE 

JULY 

576 

564  I 

728 

740 

676 

634 

524 

366 

484 

526 

JULY 

AUG. 

428 

446 

650 

702 

638 

572 

438 

376 

458 

496 

AUG. 

SEPT. 

392 

498 

516 

462 

54o 

ao8 

446 - 

364 

490 

- 438 

SEPT. 

OCT. 

_ 

356 

-  516 

340 

532 

354 

280 

296 

306 

352 

OCT. 

NOV. 

276 

290 

328 

300 

288 

294 

252 

340 

312 

234 

NOV. 

DEC. 

222 

_ 446  . 

_ m _ 

388 

344 

L55Q] 

_ zm _ 

348 

356 

_ .316 _ 

DEC. 

A|3- VALUES 


JAN. 

436 

57U 

516 

504 

.  516 

428 

378 

366 

426 

JAN. 

FEB. 

426 

580 

514 

496 

508 

412 

366 

372 

426 

FEB. 

MAR. 

jao , 

586 

-.-  500 

484 

490 

402 

360 

382 

426 

MAR. 

APR. 

430 

CM 

00 

tr\ 

486 

490 

476 

392 

348 

376 

414 

APR, 

MAY 

428 

582 

474 

486 

458 

384 

354 

378 

4lO 

MAY 

JUNE 

442 

526 

476 

488 

462 

384 

362 

378 

4lO 

JUNE 

JULY 

388 

456  “ 

592 

T70 

488 

— 552 

384 

362 

380 

406 

JULY 

AUG. 

402 

466 

590 

466 

494 

464 

382 

354 

378 

4oo 

AUG. 

SEPT. 

226 , 

474 

584 

462 

VJ1 

0 

0 

45S 

382 

356 

388 

402 

SEPT. 

OCT. 

400 

530 

590 

486 

5l4 

454 

390' 

358 

4o4 

402 

OCT. 

NOV. 

432 

562 

554 

494 

510 

444 

388 

354 

412 

398 

NOV. 

DEC. 

438 

566 

530 

sgu 

510 

”486  “ 

388 

328 

420 

378 

DEC. 

I -VALUES 


JAR. 

0.674 

0-697 

0-736 

0.813 

0,729 

1.925 

0.905 

0.814 

0.836 

JAN. 

FEB. 

0.887 

0.748 

0.732 

0.685 

0.953 

0.942 

0.880 

0.629 

0.690 

FEB. 

MAR. 

0.656 

0.802 

O.728 

0.921 

O.878 

0.955 

1.056 

0.911 

0.836 

MAR. 

APR. 

1.019 

1.736 

1.111 

i-3<*3 

1.076 

1.041 

i74i4 

1.106 

1.328 

APR. 

MAY 

2.028 

1.419 

1.190 

..  1-305 

1.306 

1.010 

1.096 

1.212 

1.283 

MAY 

JUNE 

1.567 

1.122 

1.156 

1.381 

1.292 

1.055 

1.153 

1.327 

JUNE 

JULY 

1.484 

1.237 

1.230 

1.374 

1.385 

1.372 

1.365 

1.011 

1.274 

1.296 

JULY 

AUG. 

1.065 

0.957 

1.102 

1.506 

1.291 

 1-233  - 

1.147 

1.062 

1.212 

1.240 

AUG. 

SEPT. 

0.990 

1.051 

0.884 

1.000 

1.080 

o.82i  

I.168 

1.022 

1.263 

1.090 

SEPT. 

OCT. 

0.850 

0.672 

0.873 

0.700 

1*035.  - 

0.780 

0.718 

0.827 

0.757 

0.876 

OCT. 

NOV. 

0.639 

0.5l6 

0.682 

0.607 

0-565 

OT662 

0.649 

0.960 

0.757 

0.588 

NOV. 

DEC. 

0.507 

0.788 

0.921 

0-783 

0.674 

‘  TT.261T 

0.747 

0.972 

0.848 

0.836 

DEC. 

[  ]t  "cans  data  were  missing  and  for  purposes  of  continuity  have  been 
interpolated  by  the  authors. 
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Table  A~l.  — Slough—  Continued 

YEAR 

I  1945  ‘l  1946  1  1947  I  1948  I  1949  I  1950  I  1951  I  1952  ]  1933  1 


A -VALUES 


JAN. 

mo 

mo 

44o 

540 

55°  

— 

— mo — 

m 

_ 

FES. 

280 

— m — 

490 

400 

44o 

47O.. 

486 

590 

360 

MAR. 

320 

420" 7 

510 

456 

410 

390 

310 

350 

395 

APR. 

353 

mo 

460 

530 

455 

440 

465 

370 

430 

MAY 

370 

1+40 

580 

560 

435 

446 

500 

4io 

400 

JUNE 

390 

480 

570 

480 

44o 

1*15 

520 

420 

JULY 

370 

510 

550 

520 

470 

420 

520 

4i5 

AUG. 

360 

510 

590 

520 

490 

460 

515 

440 

SEPT. 

350 

460 

560 

445 

415  . 

48o 

480 

 390 

OCT. 

330 

350 

4oo 

3S5 

410 

Too 

340 

 35j 

NOV. 

310 

360 

350 

380 

390 

415 

410 

346 

DEC. 

_ 220 _ 

360 

390 _ 

460 

410 

530 

_ m _ 

330 

ai3-  VALUES 


JAN. 

364 

400 

472 

485> 

455 

425 

458 

451 

FEB. 

354 

412 

480 

483~ 

452 

424 

466 

445 

MAR. 

342 

420 

482 

472 

444 

423 

467 

* 425 

APR. 

334 

420 

479 

458 

442 

422 

T57 

426 

MAY 

330 

422 

479 

457 

S42-- 

423 

457 

425 

JUNE 

342 

426 

480 

465 

444 

533 

463 

420“ 

JULY 

348 

430 

494 

478 

437 

437 

472 

420 

AUG. 

358 

436 

492 

470 

431 

TSS 

”475 

398 

SEPT. 

370 

hkT~ 

488 

471 

427 

+3+ 

465 

391 

CCT. 

376 

446 

490 

¥72  1 

4^  1 

44o 

469 

397 

NOV. 

382 

458 

49T  “ 

465 

1+28 

444 

465 

4oo 

DEC. 

392 

468 

489 

_ _ 

] 427 

450 

m 

400 

I- VALUES 


I- VALUES 


JAN. 

0.990  

FEB. 

0.886 

MAR. 

O.881 

APR. 

1.138 

MAY 

1.208 

JUNE 

1.180 

JULY 

1.211 

AUG. 

1.147 

SEPT. 

1.035 

OCT. 

07826 

NOV. 

0.746 

DEC. 

67884 

JAN. 

1.126 

1.025 

0.932 

 1-134 

1.209 

67859 

1.004 

1.430 

FEB. 

0.791 

1.068 

1.021 

67828 

0.973 

I7i68~ 

1.030 

1.101 

MAR. 

0.936 

1.000 

1.658 

0.932 

0.923 

0.922 

0.664 

0.805 

APR. 

1.078 

1.000 

0.960 

1.157"  . 

1.029 

17643 

1.018 

0.869 

MAY 

1.121 

1.043 

1.211 

1.225 

0,984 

l.o4o 

1.094 

0.962 

JUNE 

1.146 

1.127 

1.188 

1.032 

..  0.991  . 

0.958 

1.123 

1.012 

JULY 

1.063 

1.186 

1.113 

1.088 

1.076 

0.961 

1.102 

0.988 

AUG. 

1.006 

I.I70 

1.199 

1.106 

1.137 

1.031 

1.084 

1.106 

SEPT. 

“  67646 

i.o4i 

1.148 

0.945 

0.972 

1.106 

1.032 

0.997 

OCT. 

6 !  878 

0.785 

67816 

0.816 

0.956 

0.909 

0.725 

“67854 

NOV. 

0.812 

0.786  "" 

0.706 

0^17 

6.911 

0.935 

0.882 

0.850 

DEC. 

0.944 

0.769 

0.798 

1.011 

0.960 

1.178 

1.078 

67825 
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Table  A— II.  —  Freiburg 
YEAR 


1  19L8  |  1949  I  1950  I  1951  I  1932  |  1953  1  1954  I  1955  1  19^6  |  1957  I 


A  -  VALUES 


JAN. 

- U2B - 

331 

286 

380 

2S0 

290 

190 

430 

516 

JAN. 

FEB. 

360 

555 

222 

390 

260 

250 

250 

390 

420 

FEB. 

MAN. 

380 

J36 

250 

410 

260 

280 

280 

400 

420 

MAN. 

APR. 

m 

450 

380 

270 

286 

0 

00 

on 

470 

600 

APR. 

MAY 

j+5Q 

500 

400 

340 

410 

360 

420 

600 

MAY 

JUNE 

592 

450 

516 

310 

320 

280 - 

570 

690 

JUNE 

JULY 

5^2 

480 

jflO 

300 

320 

270 

640 

660 

JULY 

AUG. 

519 

570 

440 

300 

310 

446 

56o~ 

556 

AUG. 

SEPT. 

481 

531 

Too 

280 

300 

340 

520 

656 

SEPT. 

OCT. 

322 

280 

260  . 

^ 20Q  , 

310 

360 

490 

540 

OCT, 

NOV. 

330 

336 

360 

200 

240 

300 

360 

470 

510 

NOV. 

DEC. 

365 

361 

370 

290 

200 

310 

466 

~w 

670 

DEC. 

ai3-  values 


JAN, 

467 

1+25J 

rpr 

314 

282 

297 

422 

518— 

JAN. 

FEB. 

468 

[422] 

[181. 

30 1 

283 

306 

432 

510 

FEB. 

MAR. 

jt£] 

[421] 

[779; 

288 

...  283 

308 

438 

529 

MAR. 

APR. 

444 

421] 

771 

273 

285 

312 

450 

531 

APR. 

NAY 

451 

447 

[520] 

[369. 

271 

293 

317 

- - 558- 

532 

MAY 

JUNE 

423 

450 

>19! 

'  w 

271 

298 

324 

562 

548 

JUNE 

JULY 

451 

444 

[418’ 

365 

271  

298 

333 

463 

566 

IULY 

AUG. 

455 

1+1+1+ 

LOU 

326 

269 

294 

P9 

562 

589 

AUG. 

SEPT. 

J+53 

429 

T09] 

346 

.  .  271  . 

297 

361 

464 

609 

SEPT. 

OCT. 

^59 

428 

'406 

335 

272 

304 

375 

480 

63? 

OCT. 

NOV. 

1+68 

IT261 

898" 

332 

283 

310 

378 

490 

636 

NOV. 

DEC. 

468 

T+251 

322 

Ibi 

301 

395 

ao 

651 

DEC. 

I- VALUES 


JAN. 

0.709 

To. 659' 

0.980" 

0.892 

1.028 

0.640 

1.019 

0.791 

JAN. 

FEB, 

1.036 

[0.924" 

[-1.0081 

0.864 

0.883 

0.817 

0.903 

0.824 

FEB. 

MAR. 

0.729 

[0.591+] 

[1.082] 

0.903 

0.989 

0.909 

0.913 

0.794 

MAR. 

APR. 

1.014 

T.024; 

0.989 

0.982 

1.218 

1.044 

1.130 

APR, 

MAY 

0.998 

1.119 

[1.186, 

1,255 

1.399 

1.136 

0.917 

1.128 

MAY 

JUNE 

1.086 

1.000 

'1.111 

1.144 

1.074 

0.864 

1 . 234 

1.259 

JUNE 

JULY 

1.224 

1.081 

1.288 

1.107 

1.074 

~o78ii 

1.382 

1.166 

JULY 

AUG. 

1.141 

1.284 

1.236 

1.115 

1.054 

1.261 

0.866 

0.917 

AUG. 

SEPT. 

1.062 

—  1‘.005  j 

. .  1-156 

 U033 

1.010 

0.942 

1.121 

1.067 

SEPT. 

OCT. 

0.702 

7o.65tT 

0.776 

0.735 

1.020 

0.960 

1.021 

0.852 

OCT. 

NOV. 

0.718 

.0.845] 

0.602 

0.848 

0.968 

0.952 

0.959 

0!  802 

NOV. 

DEC. 

0-771 

.0.868] 

0.901 

0.712 

1  1.030 

1.013 

o.8o5 

1.029 

DEC, 
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Table  A-II.  —  Freiburg— Continued 
YEAR 

T"  1958  |  1959  I  I960  |  1961  I  1962"  I 


I  -  VALUES 


JAN. 

0.876 

FEB. 

0.901 

MAR. 

0.881 

APR. 

1.086 

MAY 

1.142 

JUNE 

1.140 

JULY 

1.160 

AUG. 

1.104  ' 

SEPT. 

1.063 

OCT. 

o7B38 

NOV. 

o78o6 

DEC. 

0.876 
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Table  A— III.  —  Dakar 
YEAR 

r  i  1931  1  1932  1  1953  1  195^  1  1955  1 


A -VALUES 


JAN. 

310 

220 

210 

146 

FEB. 

i4o 

170 

260 

170 

MAR. 

120 

310 

290 

170 

APR. 

400 

440 

330 

190 

MAY 

370 

552 

350 

'  450 - 

170 

JUNE 

450 

210 

100 

310 

200 

JULY 

310 

100 

360 

350 

180 

AUG. 

34o 

120 

 330 

350 

210 

SEPT. 

186 

456 

390 

270 

310 

OCT. 

136 

4  80 

370 

200 

380 

NOV. 

210 

236 

230 

230 

0 

on 

on 

DEC. 

270 

210 

_  ,  2;o _ 

260 

220 

A|3-  VALUES 


JAN. 

256 

282 

315 

222 

FEB, 

241 

293 

314 

?n 

MAR. 

249 

315 

309 

208 

APR. 

272 

303 

295 

216 

MAY 

2I2 

289 

284T 

226 

JUNE 

279 

284 

286 

225 

JULY 

216 

285 

281 

AUG, 

265 

287 

278 

SEPT. 

 278 

296 

271 

OCT, 

302 

298 

263 

NOV. 

283 

222 

298 

251 

DEC. 

276 

271 

2.95 

232 

I -VALUES 


JAN. 

1.211 

.780 

.667 

.631 

FEB. 

.581 

1 580 

~  ^828 - 

^Bo8 

MAR. 

.482 

.  .987 

.939 

.817 

APR. 

OTu  . 

l74f9 

1.119 

OTo - 

MAY 

1.153 

1.311 

1.585 

.752 

JUNE 

.753 

.352 

1.084 

.889 

JULY 

^688 

1.268 

1.246 

AUG. 

.453 

1.150 

1.259 

SEPT. 

1.583 

1.318 

.996 

OCT. 

1.589 

1.242 

.760 

NOV. 

.729 

.936 

.772 

.916 

DEC. 

^28 

775 

.712 

1.121 

I- VALUES 


JAN. 

.822 

FEB. 

.699 

MAR. 

.806 

APR. 

1.225 

MAY 

1.326 

JUNE 

.770 

JULY 

-1 .067 

AUG. 

t954 

SEPT. 

1.299 

OCT. 

\  1.197 

NOV. 

.MS 

DEC. 

00 
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Table  A-IV.  —  Singapore 


YEAR 

|  1949  |  1950  1  1951  |  1952  I  1953  1  1954 j 


A- VALUES 


JAN. 

^8o 

'lio] 

315 

120 

280 

FEB. 

405 

[360] 

315 

330 

305 

MAR. 

470 

390 

120 

370 

120 

APR. 

520 

405 

110 

-150 

270 

MAY 

500 

[375] 

110 

100 

JUNE 

125  

150  

 r^i 

270 

JULY 

330. 

310 

300 

270 

AUG. 

400 

160 

3U0 

105 

SEPT. 

515 

380 

400 

380 

350 

OCT. 

3k5 

175 

425 

400 

365 

NOY. 

505  ~ 

350 

380 

366““ 

295 

DEC. 

410 

100 

110 

_ m _ 

JML 

A|3- VALUES 


JAN. 

358  . 

117 

112 

— 

FEB. 

160 

-.119 _ 

333 

MAR. 

432 

3jft 

333 

APR. 

42d 

384 

341 

332 

MAY 

412 

384 

336 

324 

JUNE 

396 

161 

332 

119 

JULY 

..390 

362 

333 

315 

AUG. 

188 

361  

S34 

314 

SEPT. 

357 

358 

338^ 

313 

OCT. 

182 

353 

34o 

305 

NOV. 

371 

34b 

338 

DEC. 

_ m _ 

341 

_ m _ 

I- VALUES 


JAN. 

.922 

.915 

.964 

FEB. 

1.000 

-929 _ 

>991 

MAR. 

1.076 

1.081 

.918 

1.111 

APR. 

1.221 

1.113 

.968 

1.054 

MAY 

1.214 

1.030 

.923 

7926 

JUNE 

7821 

.970 

«919 

.846 

JULY 

.846 

.856 

.901 

.857 

AUG. 

I.031 

.997 

1.018 

•  971 

SEPT. 

.982 

1.117 

1.124 

1.118 

OCT. 

.982 

1.204 

1.176 

1.197 

NOV. 

•  9^3 

17698 

1.065 

DEC. 

.816 

.909 

.985 

I -VALUES 


JAN. 

.940 

FEB. 

.973 

MAR. 

1.052 

APR. 

1.089 

MAY 

1.7)23 

JUNE 

7889 

JULY 

78^5 

AUG. 

1.004 

SEPT, 

I.O85 

OCT. 

1.140 

NOV. 

1.035 

DEC. 

.910 
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